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Abstract 

The HyperCP collaboration has recently reported the observation of three events for the decay S + — * 
pfi + fi~ with an invariant mass m^+^- for the muon-antimuon pair of ~214MeV. They suggest that 
a new particle state X may be needed to explain the observed m„+„- distribution. Motivated by 
this result, we study the properties of such a hypothetical particle. We first use K + — > 7r + /i + /x~ 
data to conclude that X cannot be a scalar or vector particle. We then collect existing constraints 
on a pseudoscalar or axial-vector X and find that these possibilities are still allowed as explanations 
for the HyperCP data. Finally we assume that the HyperCP data is indeed explained by a new 
pseudoscalar or axial-vector particle and use this to predict enhanced rates for Kl — > ttttX — ► 7T7r/i + /x _ 
and f2~ — > 3~X — > 
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I. INTRODUCTION 



Three events for the decay mode S + — > pfi + pT with an invariant mass of 214.3 ± 0.5 MeV 
for the muon-antimuon pair have been recently observed by the HyperCP collaboration (lj . The 
branching ratio is obtained to be [8.6l|;4(stat) ± 5.5(syst)] x 1CT 8 [1]. The central value is 
considerably larger than the short-distance contribution in the standard model 0. When long- 
distance contributions are properly included, it is possible to account for the total branching 
ratio 0, El- However, the clustering of the events [whose contribution to the branching ratio 
is (3.lli'g ±1.5) x 10 -8 ] around 214 MeV cannot be explained. If this result stands future 
experimental scrutiny, it is most likely to be due to a particle state X having a mass of 214 MeV. 
In this paper we study the properties of such a particle assuming its existence. 

The mass 214 MeV of this hypothetical particle is close to, but higher than, the sum of the 
masses of two muons. It is tempting to identify it as a muonium bound-state. However, the 
S-wave bound-state has a mass below the sum of the two muon masses. Therefore, the state 
X cannot be an S-wave muonium state. Radial excitations can yield larger masses, but it is 
unlikely that the electromagnetic interaction bounding the muon and antimuon together can 
raise the mass by the 3 MeV needed. The X particle, if exists, is likely a new state beyond the 
standard model (SM). There are theories where such light states naturally exist, for example, the 
super-partner of the goldstino particle in spontaneously local super-symmetry breaking theories 
as discussed in Ref. |4J. These particles, pseudoscalar and scalar ones, can have masses lower 
than a few GeV or even in the MeV range. 

In our study we will not attempt to construct models which predict such particles. Instead, 
we will assume the existence of the new particles and study the implications from the HyperCP 
data. To be consistent with observations, we follow HyperCP and assume that the hypothetical 
particles have small widths, are short-lived (they decay inside the detector), and do not interact 
strongly [lj]. 

II. EFFECTIVE INTERACTIONS 

In our study, we will try to be as model independent as possible by parameterizing the 
interactions of this new particle with known particles. In the HyperCP hypothesis, this particle 
is produced in the decay of £ to p and subsequently decays into a muon-antimuon pair. At the 
quark level, the particle X must then couple to ds (and of course to as well). A priori, 

the state X can be a scalar, pseudoscalar, vector, axial-vector, or even a tensor particle. We will 
consider four possibilities: scalar (X$), pseudoscalar (Xp), vector (.XV), and axial vector (Xa)- 

Assuming that the hypothetical new particles have definite parity, do not carry electric or 
color charge, and are their own anti-particles, we can write their couplings to ds and as 

£s = (-Qsgds + H.c.)X 5 + g Sfl finX s , 

£p = {-Wpqdl^s + H.c.)X P + ig Pl ^flj 5 p,Xp , 

£v = {-9v q dl^s + H.c.)X£ + gvppnpliXv » 

C-a = {gAqdltilbS + K.c.)X% + # A/i /i7 M 75/iX^ . (1) 

If the particle does not have a definite parity, our results should be interpreted as applying to 
the parity-even or -odd coupling as appropriate. 
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A condition that the couplings in the above equations must satisfy is that they must be able 
to produce the observed branching ratio in S + — > pfi + fi~. To carry out such a fit, one must 
know how the X couples to the hadron states S + and p from the above quark-level couplings. 
To this end, we employ chiral perturbation theory to obtain the couplings. Our task is simplified 
by the assumption that the hypothetical particles do not interact strongly as they can then be 
readily identified with the scalar, pseudoscalar, vector, and axial-vector external sources in the 
standard-model Lagrangians. With the flavor properties assumed in Eq. (fl|). the appropriate 
Lagrangians are then 



£sb v = b D (B {gh 8 ? + £hgt,B}) + b F (B [tfhgtf + thg£, B] > 

+ b Q (h s (St + £)) (BB) + \fB, (h s (St + S)> + H.c. , 



(2a) 



£ 



PBip 



ib D (B {Ch p e - £h P Z,B}) + ibp (B [?h P ? - Shpt, B] ) 
+ ib (h P (tf-Z))(BB) + ~fB (h P (tf-Z)) + H.c. 



(2b) 



£vb v = \{Bl,[B,eK^ + ^\ 

+ \d (b %15 {^h^ - b}) + 1 -f (b %15 [fhitf - e^e f , B] 

+ \C[% i^Ki - £h!$) B + B {^Ki - ih^) T„ ] 
- ^/ 2 <^(^SSt-St^S)> + H.c. , 



(2c) 



£ 



ABip 



B, tfhft - thtf 



+ \d (b%^ {Wa* + #itf f . B }) + \ F ( S W* [fhft + frtf, B 



+ \ C [% + ^) B + B {tfhft + ] 

- H.c. , 



(2d) 



where we have shown only the terms relevant for this paper, and used the notation (h Y ) kl = 
X Y9 Yq Ki for Y = S,P, and (h Y ) kl = X Y g Yq h kl for Y = V, A, with h kl = (T 6 + iT 7 ) kl = 
S k2 S 3l . The notation and parameter values that we employ here are explained in Appendix 1X1 
With the above effective Lagrangians, we can obtain constraints on the couplings g Yq from other 
low-energy processes. 

III. RULING OUT THE SCALAR AND VECTOR AS CANDIDATE PARTICLES 



With the assumption that the new particles are short-lived and narrow, their contribution to 
the branching ratio of S + — ► pfi + fi~ is given by 23(S — > pX)B(X — > /i + /i~). Using the effective 
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Lagrangians in Eq. we find the matrix elements for S + — > pX to be 

JU(E+ - P X P ) = g Pq B (D - F) m 2 S + mp 2 p 75 S + , 

y — nip 

.M(E+^pX y ) = -g Vq pfV% , 

M(Z+ - pX A ) = -g Aq (D - F)^ 75 S+ e ; . (3) 

These expressions follow from a kaon-pole diagram for the pseudoscalar, and from a direct vertex 
from Eq. (J2J) for the rest. For the branching ratios, it then follows that 

= 9.0 x W> X2 \g Sq \ 2 B{X s ^^iT) , 
= 3.7xlO n \g Pq \ 2 B(X P ^fi + fx-) , 

= 7.0 x 10 1:L \g Vq \ 2 B(X V fi + {i~) , 
= 7.0xl0 11 \g A fB(X A ^ f i+^) . (4) 





^pX s - 




> 


- 


^pX P - 




> 


- 


^pX v - 


■*p(f 


> 


- 


^pX A - 




> 



For the scalar and vector particles, there are severe constraints from — > 7r ± /i + /i~. The 
branching ratio of K ± — > n ± fi + fi~ has been measured to be B = (8.1 ± 1.4) x 10~ 8 |5[. The X- 
particle contribution to these decays can again be factorized as BiK^ — > n ± X)B(X — ■> ^). 
Using the effective Lagrangians in Eq. (J2J), we have the matrix elements for K — > ir X 

M^^^Xs) = g Sq B , 

MiK* - 7r ± X y ) = g Vq (p K + Pv ) ■ e* . (5) 

We have assumed CP conservation for simplicity, and so taken the couplings g, s V s q to be real. 

The decay modes K — > 7r ± /i + /i~ are long-distance dominated in the SM [6] and the measured 
spectra agree reasonably well with the predictions 0, 0] . In particular, there is no apparent bump 
in the m^,,- = 214 MeV region jjj. In view of this, we require that any contribution from the 
hypothetical new particles to these rates be below the experimental error, that is 

BiK* -> vr ± /i + / u-) x < 1.4 x 1(T 8 . (6) 

This leads to the constraints 

\g Sq \*B{X s ^ii+vr) < 6.5xl(T 24 , \g Vq \ 2 B(X V -> < 4.3 x 1(T 23 . (7) 

Combining these limits with Eq. (|4*j). we find 

£(E+ -»■ pX s -> p/i + /0 < 6 x 10 -11 , -> pX y -> p/i + /0 < 3 x KT 11 . (8) 

These results indicate that — > 7r data rule out both a scalar particle and a vector 

particle as explanations for the HyperCP result. Notice that this conclusion still holds if we relax 
Eq. © and allow the new contribution to be as large as the full experimental rate. 

The decays K ± — > n ± Xp i A_ are not allowed, as we have assumed X PA to have no parity- 
odd couplings. Therefore, there are no constraints from K — > TTfi + fi~ on the pseudoscalar and 
axial- vector couplings of the hypothetical particles to quarks. 
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IV. SOME CONSTRAINTS ON PSEUDOSCALAR AND AXIAL- VECTOR COU- 
PLINGS 



We now consider other possible constraints on the couplings involving the pseudoscalar and 
axial- vector particles. We begin by ignoring CP violation so that g Pq and g Aq are real. A strong 
constraint on flavor-changing neutral currents (FCNC) of this type comes from K°-K° mixing. 
The mixing parameter M 12 = Ai(K° — > X — > K°)/2mK from an intermediate X-state in 
K° ^ X -> K° is given by 



M{K° ^X P ^K L 



2B* f*g% q 

2 2 

m K — m p 



M{K o ^ Xa ^ K°) = ^ . (9) 



m A 



The measured value of AMk l -k s — 3.483 x 10~ 12 MeV |5j can be accommodated in the SM, 
but its calculation suffers from hadronic uncertainties due to long-distance contributions. To 
be conservative, we will thus require that any new physics contribution be smaller than the 
experimental value, namely 

(AM Kl _ Ks ) x = 2 (ReM 12 ) x < 3.483 x 10" 12 MeV . (10) 



With matrix elements from Eqs. (|2b|) and (|2dj). but using f K ~ 1.23/, instead of /, in Eq. (J9J) 
for the kaon decay constant, this results in 

g 2 Pq < 3.3 x 10- 15 , 

g\ q < 1.3 x 10- 14 . (11) 
When we substitute these bounds into Eq. (JH), we find 

£>(£ + — ► pXp — > pfi + fj, 



B(X P -> H+/JL-) 
S(S + — > pX A — > P^ + ^~ 
B(X A ^ fi+fj,-) 



< 1.2 x 10 -3 



< 9.1 x 1(T 3 . (12) 



These constraints are so weak that Xp and X A are allowed candidates to explain the HyperCP 
result, provided their branching ratios into muon pairs are at least B{X P — > fi + fi~) > 2.5 x 1CT 5 
and B(X A — > > 3.4 x 10 -6 , respectively. 

If we allow for CP violation in the g Pq and g Aq couplings, the constraints from AM Kl _ Ks 
become 



-15 



|(Re^) 2 - (Im^) 2 ! < 3.3 x 10 . 

-18 



\(Reg Pq )(lmg Pq )\ < 3.2x10" . 
|(Re^ Ag ) 2 - (Im^) 2 | < 1.3 x 10" 14 , 

|(Re^)(Im^ ? )| < 1.2 x 10- 17 . (13) 

The two additional constraints arise from the new-particle contribution to the parameter = 
IH1M12/ (v2 AMk l -k s )- This parameter can be calculated more reliably than AMk l -k s in the 



5 



standard model and is in good agreement with the result \ex\ = 2.284 x 10~ 3 . In view of this, 
we required the new-physics contribution to be less than 30% of the experimental value, which 
is about the size of the theoretical uncertainty in the SM calculation. 

New pseudoscalar and axial-vector particles also contribute to the rare decay Kl — > 
via the pole diagram Kl — > X — > fx + fi~ . From Eqs.(JH) and (0), we obtain 

-2iB fg Pq 



M(K L 
M(K L 

These matrix elements imply 



B[K L 
B(K L 



X A 



x f 



X, 



fx fx 

fX + fX~ 



fx fx 



fl^fX 



2 2 9pnV>l5V 

m K — mp * 



(14) 



5.6 x 10 18 GeV" 1 g\ T(X P 
1.2 x 10 18 GeV^ 1 g% T(X A 



fx fl 

fX + fX~ 



(15) 



If we allow for CP violation it is possible to obtain additional, weaker, constraints from consider- 
ing the mode Kg — > fx + fxT . To be useful, the equations above must be combined with additional 
information on the couplings of the hypothetical new particles to muons. Partial information 
can be obtained from considering their contribution to the anomalous magnetic moment of the 
muon, a M . 

At one-loop level, the contributions of the new pseudoscalar and axial-vector to are given 
respectively by 



a^A) 



\9p. 



\ 2 m 2 



\9al 



-2.28 x 10 



-3 



9pl 



47r 2 m A 



f A {m 2 Jm 2 A ) = -8.97 x 10" 3 \g A . 



(16) 



Here 



fp(r) 
f A (r) 



dx 



dx 



x 



1 — x + rx 2 ' 

4(x — l)x + x 2 (l — x) — 2rx 3 



(17) 



1 — x + rx 2 

At present there is a discrepancy of 2.4cr between the SM prediction and data [10], Aa^ = 
a^exp) - a^(SM) = (23.9 ± 10) x 10~ 10 with a^(exp) = (11659208 ± 6) x 10" 10 . We note that 
the new contributions reduce the value of a M , making the comparison with experiment worse. In 
view of this, we place a conservative constraint on g\ by requiring that the new contribution 
to a M not exceed the experimental error. This results in 



\9Pft 
\9 A ^, 



2 < 2.6 x 10~ 7 
2 < 6.7 x 10" 8 



r(x P 
r(x A 



fTfX 
fX + fX~ 



< 3.7 x 10~ 10 GeV 



< 5.2 x 10~ 12 GeV . 



(18) 



Combining the constraints in Eqs. (fTT|) and (fTSj) . we obtain from Eq. 



B{K L 
B(K L 



X P 
X A 



fx + fi-) < 6.8 x 10" 
fx + fi-) < 8.1 x 10" 



(19) 
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The measured branching ratio for this mode, B(Kl — > ^ + ^) = (6.87 ± 0.12) x 10~ 9 , is 
almost completely saturated by the two-photon intermediate state, the absorptive part of this 
contribution being B(K L — > 77 — > /i + /U~) abs = (6.63 ± 0.07) x 10~ 9 (this is referred to as the 
unitarity bound). This leaves little room for a direct new-physics contribution. Here we assume 
that a possible new-physics contribution is at most equal to the difference between the measured 
rate and the unitarity bound plus one standard deviation, 

B(K L -> fi + fi~)x < 3.6 x 10" 10 . (20) 

Using this as a constraint improves the bounds of Eqs. (fTTJl and (fTHj) : 

\g Pq \ 2 T(X P -> fi + fT) < 6.4 x 10~ 29 GeV , 

\g A /T(X A -> fi+fT) < 3.0 x 10" 28 GeV . (21) 

If T(Xp t A — > H + H~) are allowed to saturate the bounds of Eq. (|TH|) . the above equations imply 
that 

\g Pg \ 2 < 1.7 xlO" 19 and therefore p// + /0 < 6.3 x 10~ 8 B(X P -> , 

|c/ A / < 5.8xl0~ 17 and therefore £(£+ -> p/x + /0 < 4.0 x 10~ 5 i3(X A -> // + /0 . (22) 

This in turn means that both the pseudoscalar and axial- vector particles remain viable candidates 
to explain the HyperCP data after combining the existing bounds from AMk l -k S j a M' anc ^ 
K L — > fi + fi~. In the case of the pseudo-scalar, these combined bounds require that it decay 
almost exclusively into a yU + /i~ pair. 



V. PREDICTIONS 



We now turn the argument around and assume that the HyperCP data is indeed explained 
by the hypothetical new pseudoscalar or axial-vector particle. This implies that 

\g Pq \ 2 B(X P ^ f i+^) = (8.4+^±4.l)xl0- 20 , 

\g Aq \ 2 B(X A ^^^) = (4.4^;* ±2.1) x lO" 20 . (23) 

These can then be used to predict the contributions of the new particles to other decay modes 
such as K — » TT7rXp yA — > 7T7r/i + /i~ and Q~ — > E~X PtA — > S~/i + //~. 

We first consider K — > ititX Pj a — > irir fi + . Employing the Lagrangians in Eqs. and (53) , 
we derive 1 

.M(X U -> TT + TT-Xp) 
7W(i?° ^7T 7T Xp) 



^Pg + ™* 



Ci -+--\~ _ t,i, -''/'</ l"'/v '"rr _m X7r- ~ m *+7r- 



2 2 



£oflp g {™? K -m\- m 2 07r0 ) 

2V27 



2 2 
m A" — m P 



W -> 7T + 7T Xp = -^p- ££L , 24 

2/ m^- — mp 



1 We note that each of the K° — > 7T7rXp amplitudes receives contributions from both contact and kaon-pole 
diagrams. The pole terms seem to be missing in Ref. 
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M(K° - Ti+71- X A ) = -il^Lp^ ■ e * , 
M(K°^n\°X A ) = -i%p K -e\ 

M (K + - 7r + n°X A ) = -i^ip^o - Pw+ ) • e* , (25) 

where m?- = (p^ + £>j) 2 - 

Adding the errors in Eq. (|23|) in quadrature, we obtain the predictions 

B(K L -> TT+TTXp -> 7T + 7r> + /0 = (l.8±^) X 1(T 9 , 

B(Kl -> vrVXp -> ttVV/O = (8.3j£jj) x 10" 9 , (26) 

£(X L -> 7T + 7T-X A -> Tr+TT/X+pT) = (7.3±|:f) X 10" 12 , 

B (K L - vrVX^ - 7T°7r W) = (l.0+°;i) x lO' 10 . (27) 

Notice that these decay modes are highly suppressed by phase space, but that at the lCT 8 -iCr 9 
level they are comparable to existing limits on other rare Kl decay modes. The rates for the 
K + decay modes are quite sensitive to isospin-breaking effects, but we find them to be at most 
at the 10 -12 level, much less promising than the modes. 

We now consider the modes Q~ — > E~Xp }A — > E~ . For the pseudoscalar particle, a 
kaon-pole diagram with vertices from Eqs. (jAlj) and (|2b|) leads to 



M (0- - ETXp) = J B ' C g Pq u s q^ Q , (28) 

K P 

and for the axial-vector particle a direct vertex from Eq. (|2d[) gives 

M (O- - E~X A ) = -iC g Aq « 3 < e; . (29) 

The resulting branching ratios for (Q~ — > S~X — > E~fi + fi~) are 

1 Id I 3 B 2 C 2 \a I 2 
B(QT -> E-Xp -> E-fi+fj,-) = ^ ^^(E = + m = )i3(X P ^/i + /i") 

= 2.4 x 10 13 |s Pg | 2 £(P^ + pT) , (30) 
B(fi- -> S-X A -> E-fi+fx-) = ^- c 2 \g Aq \ 2 (E B + ma) f 3 + 4") - /iV) 



T c _ 12vrm Q 



= 1.6 x 10 13 |^ 9 | 2 S(X A -> //+/0 . (31) 

Consequently, the HyperCP data implies 

B(fi- -> E-Xp -► H>+/0 = (2.0±i;l ± 1.0) x 10~ 6 , 

-> 2~X A -> E~/j, + h-) = (0.73+°;^ ± 0.35) x 1(T 6 . (32) 

These numbers represent a substantial enhancement over the existing standard-model prediction 
B(p,~ -> ETfi+fj,-) = 6.6 x lO- 



ll 
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VI. SUMMARY 



We have studied the hypothesis that a new particle of mass 214.3 ± 0.5 MeV is responsible 
for the invariant-mass rriu+u- distribution observed by HyperCP in S + — ► pfi + fi~ . We find that 
existing data on K + — > 7r + /i + //~ rule out a scalar particle and a vector particle as possible 
explanations. We explore all the existing constraints on pseudoscalar and axial-vector particles, 
and conclude that these possibilities are still allowed. If either one of them is indeed responsible 
for the HyperCP data, we predict enhanced rates for Kl — > ttttX — ► 7T7Tfi + fi~ and Q~ — ► 
S-X -> E-/j,+/j,-. 

Note added After the completion of our paper, the work of Deshpande, Eilam and Jiang 
appeared. They reach similar conclusions to ours. 
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APPENDIX A: DERIVATION OF EFFECTIVE LAGRANGIANS 

The chiral Lagrangian that describes the interactions of the lowest-lying mesons and baryons 
is written down in terms of the lightest meson-octet, baryon-octet, and baryon-decuplet fields [13, 
0- Q|- The meson and baryon octets are collected into 3x3 matrices <f and B, respectively, 
and the decuplet fields are represented by the Rarita-Schwinger tensor T^ bc , which is completely 
symmetric in its SU(3) indices (a,b,c). The octet mesons enter through the exponential X = 
£ 2 = exp(i<f/f), where / = f n = 92.4 MeV is the pion-decay constant. 

We write the strong chiral Lagrangian at leading order in the derivative and m s expansions 

as 

£ s = (Bir(d,B+[V^B]))-m (BB) + D(BY%{A^B})+F(BYl 5 [A^B}) 
- T M ilpT^ + m T T% + C (T^A^B + BA^T") + H T^ 7s 

+k {B{x +- B})+ k {B[x+ - B])+ k M{BB) 

+ Wj"^" - Jb„ T " T " + if <jD " E ' + * f < x+) • (A1) 

where (■ ■ ■ ) = Tr(- ■ ■ ) in flavor-SU(3) space, m and m T are the octet-baryon and decuplet- 
baryon masses in the chiral limit, respectively, V = §(£<9 M £ f + f f ) + f(^ f ^ + £r"£t), 

a» = | (e <n f - e f <n) + § (e f n - ^) , v^ lm = d^ lm + v^ mn + v? n n mn + v^ n , 

D^E = <9"E+iPE-iEr", x+ = tfx^Hx^, with P = \\J» = v"+a^ = \\ a r£ = v^-a 1 *, 
and x = s + W containing external sources. In the absence of external sources, \ reduces to the 
mass matrix \ = 2B diag(m, m, m s ) = diag(m 2 , m 2 , 2m 2 K — m 2 ) in the isospin-symmetric limit 
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m u = m d = rh. The constants D, F, C, H, B , b DF0 , c, and c are free parameters which can 
be fixed from data. 

To extract the couplings of the new particles X from the above Lagrangian, we identify the 
external sources with 

s = 9sq x s ( T 6 + iT 7 ) + H.c. , 
V = 9p q X P (T 6 + iT 7 ) + H.c. , 
v» = g Vq X»(T 6 + iT 7 ) + H.c. , 

^ = 9 Aq K{T Q + iT 7 ) + H.c. . (A2) 

The Lagrangians in Eq. © then follow. 

Numerically, we adopt the tree-level values D = 0.80 and F = 0.46, extracted from hyperon 
semileptonic decays, as well as \C\ — 1.7, from the strong decays T — > Bip. Furthermore, using 
m + m s = 121 MeV and isospin-symmetric values of the baryon and meson masses, we have 

b D = 0.270 , b F = -0.849 , B = 2031 MeV , (A3) 

the other parameters being irrelevant to our calculations. 
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